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Abstract 

A new mass table calculated by the relativistic mean field approach with the state-dependent BCS 
method for the pairing correlation is applied for the first time to study r-process nucleosynthesis. 
The solar r-process abundance is well reproduced within a waiting-point approximation approach. 
Using an exponential fitting procedure to find the required astrophysical conditions, the influence 
of mass uncertainty is investigated. i?-process calculations using the FRDM, ETFSI-Q and HFB-13 
mass tables have been used for that purpose. It is found that the nuclear physical uncertainty can 
significantly influence the deduced astrophysical conditions for the r-process site. In addition, the 
influence of the shell closure and shape transition have been examined in detail in the r-process 
simulations. 
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I. INTRODUCTION 



It is of the utmost interest to explore the "terra incognita" of exotic nuclei, as evidenced 
by the fact that several Radioactive Ion Beam (RIB) facilities are being upgraded, under 
construction or planned to be constructed worldwide. Such investigations of the properties of 
these exotic nuclei, which may behave very differently from the nuclei around the /3-stability 
line, result in new discoveries such as the halo phenomenon - nucleons spread like 

a thin mist around the nucleus, which can significantly increase the nuclear reaction ratio. 
Stellar nucleosynthesis processes such as the r-process |3|, which is responsible for roughly 
half of the enrichment of elements heavier than iron in the universe, also require a thorough 
understanding of the properties of exotic nuclei. Key properties like masses for example, 
determine the path that the nucleosynthesis process follow in the nuclei chart. Nevertheless, 
despite many experimental efforts, present knowledge of exotic nuclei still does not include 
much of what is required for a complete understanding of r-process nucleosynthesis. After 
the first systematic introduction to the r-process jsj] half a century ago, r-process calculations 
for a long time could only rely on the phenomenological nuclear droplet mass formula 6| 
because of the lag of both experimental and theoretical development. Fortunately, in the 
last 15 year the theoretical study of nuclear properties has made tremendous progress and 
r-process calculations 3, S Q] have been carried out based on the refined droplet model 
FRDM [l^, Hartree-Fock approach like ETFSI-Q and the very recent microscopic 
rooted Hartree-Fock Boghubov (HFB) [l2,[l3,Q. 

Despite the progress in the theoretical nuclear structure physics, mass models predictions 
(which by design concentrate on different nuclear structure aspects) still show a large de- 
viation when going to very neutron-rich nuclides, even though they have achieved similar 
quality to describe known nuclides. This is specially troublesome since the astrophysical 
scenario in which an r-process may occur is a matter of debate and all astrophysical simu- 
lations dedicated to the nature of the stellar environment depend on the input from nuclear 
physics. Mass model predictions, even in models that give similar global rms error still show 
local deviations differently. 

In principle, microscopic-rooted mass models should have a more reliable extrapolation 
to the unknown regions, therefore these studies have received more and more interest as 
evidenced by the increasing number of non-relativistic HFB investigations 1^, Q, 
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15|. 
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Based on a mass-driven fitted method, the latest HFB models have achieved a similar quality 
[rms ~ 0.7 MeV) as the phenomenological FRDM mass model for known masses. More 
recently, another microscopic-rooted approach, the relativistic mean-field (RMF) theory [l^ 
has received broad attention due to its successful description of several nuclear phenomena 
during the past years (for recent reviews, refer to Refs. 
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18|). In the framework of 



the RMF theory, the nucleons interact via the exchanges of mesons and photons. The 
corresponding large scalar and vector fields, of the order of a few hundred MeV, provide 
simple and efficient descriptions of several important phenomena such as the large spin- 
orbit splitting, the density dependence of optical potential, the observation of approximate 
pseudo-spin symmetry, etc.. Moreover, the RMF theory can reproduce well the isotopic shifts 



in the Pb region 



19|, explain naturally the origin of the pseudo-spin symmetry 20|, l21| as a 



relativistic symmetry 
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24, 



25l | and spin symmetry in the anti-nucleon spectrum 26 1. 



The first RMF mass table was reported in Ref. 27|] for 2174 even-even nuclei with 8 < 
Z < 120 but without including pairing correlations. Later on, the calculation was improved 
by adop ting a constant-gap BCS method and calculated 1200 even-even nuclei with 10 < 



Z <9l 



28[, most of which are close to the /3-stability line. More recently, using the state- 
s', the first systematic study of the ground 



dependent BCS method with a 5-force [29 
state properties of over 7000 nuclei ranging from the proton drip line to the neutron drip 



line was performed 



311 ]. Comparison of this calculation with experimental data and to the 



predictions of other mass models will be presented in more detail in Sec II. 

Considering the recent development of the microscopic mass models in both the HFB 
and RMF approach, it is very interesting to examine their applicability to an r-process 
calculation. The main goals of this paper is to explore to what extent the solar r-process 
abundance can be reproduced using the new RMF mass table 3l|] and by comparing with 
other theoretical mass models, to determine the infiuence of nuclear mass uncertainty in 
r-process calculations. The paper is organized as follows. In Sec. |Tl]the global agreement of 
the new RMF mass table with the experimental data is discussed and the RMF prediction 
in the very neutron-rich range is compared with the FRDM |l^, the ETFSI-Q H] and 



151 ] mass tables. In Sec. Illlt a short introduction to a site-independent 



the latest HFB-13 

r-process approach is given. In Sec. IIVI the new mass table is applied to reproduce the 
solar r-process abundances. In addition, the result is compared to the r-process abundances 
obtained with the predictions of the FRDM, ETFSI-Q and HFB-13 mass models. Finally 
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the summary and conclusions are given in Sec. |Vl 



II. GLOBAL BEHAVIOR OF THE NEW RMF MASS TABLE 

With about 10 parameters fitted to the ground-state properties of around 10 spherical 
nuclei, the RMF approach with the TMA parameter set is found to give a satisfactory 
description for all the nuclei in the nuclear chart. The model deviation of one-neutron 
separation energy S'„, with respect to the known experimental data can be characterized by 
the rms deviation (arnT<) 



\ 



1 " 



n 



Although a relatively large rms deviation for the absolute mass value is found for the 
RMF calculation in comparison with the FRDM and HFB-13 models, the finite differences 
in binding energies such as the practical used one-neutron separation energies Sn are well 



predicted due to the cancelation of systematic error 
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32l | . The rms deviation of Sn for the 



FRDM, ETFSI-Q, HFB-13 and RMF models with respect to experimentally values 



are 



399 keV, 528 keV, 546 keV and 654 keV, respectively. Here the comparisons include nuclei 



with Z,N > 8. Comparing the predictions of the RMF model to the known values [J 
results in discrepancies between -1.4 MeV to 1 MeV, while the difference between either the 
FRDM or the HFB-13 and the experimental data is in the range of -1.3 MeV to 0.5 MeV. 
It shows that the microscopic model such as the RMF approach can almost achieve the 
same level of accuracy for known one-neutron separation energy Sn as the phenomenological 
FRDM. For each isobaric chain with mass A, the distance between the nuclide {Z, A) and 
the nuclide (Zq, A) in the /5-stability line 3^ is defined hj e = Zq — Z with 

^° " 1.98 + 0.0155^2/3 ' 

i.e., e = stands for the most stable nuclei and e > the neutron-rich nuclei. The rms 
deviation arms of Sn as a function oi e = Zq — Z for different mass models is shown in Fig. [H 
It is remarkable that almost the same order of prediction power of Sn from the neutron- 
deficient side to the neutron-rich side is achieved for all the models, even though models like 
FRDM and HFB-13 have about 10 more free parameters than the RMF model and were 
optimized for all the known masses. While the macroscopic-microscopic mass model FRDM 
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shows the best agreement with experimental values in the neutron deficient mass region, it 
gets progressively worse when moving away from the stability line towards the neutron-rich 
side. 

All theoretical models addressing nuclei far away from the /5-stability line involve a dra- 
matic extrapolation to unknown nuclei. Thus it is interesting to examine what is the dif- 
ference of Sn predicted in the different models when going towards the neutron-rich side. 
In Fig. [21, the differences between Sn in the RMF model and those in the HFB-13 mass 
models are shown as an example. In general, most of the discrepancies between the two 
microscopic models range from -1 MeV to 1 MeV across the entire nuclear chart. Further- 
more, the Sn value is consistent with that in the HFB-13 model in the range of -0.5 MeV to 
0.5 MeV when going to the unknown region of the nuclei chart. These differences indicates 
that the extrapolation can be quite different depending on the underlying physics of the 
model. Similar conclusion can be drawn also for the comparison of the RMF and FRDM 
models. Around the = 82 shell, the RMF model predicts a systematic lower Sn- These 
different Sn predictions towards the neutron drip line affect r-process calculations and thus 
the corresponding determined astrophysical conditions. 

The evolution of the nuclear structure around the shell closures A^ = 82 and 126 is 
critical in understanding the r-process abundance distribution around the A = 130 and 
A = 195 abundance peaks. At the shell closures, the one-neutron separation energy drops, 
and thus the corresponding nucleus in the r-process path cannot absorb another neutron 
without photo-disintegration. Therefore it has to "wait" for the /3-decay to proceed, and the 
path moves closer to the valley of stability where the half-lives are longer. These isotopes 
with long half-lives serve as bottlenecks of the process where abundances accumulate and 
the abundance peaks are formed. In Fig. [3l the predicated average one-neutron separation 
energy S'2n/2 around A^ = 82 in the RMF model is displayed as a function of mass number 
for isotopes from Kr to Ba together with the experimental values and the predictions of 
the FRDM, ETFSI-Q and HFB-13 models. The dominating isotopes in the r-process path 
(defined as waiting points and discussed in Sec. llVj) are also indicated for the corresponding 
mass model. Similarly, Fig. H] presents S'2„/2 distribution around A^ = 126 shell for isotopes 
from Ce to Pt. In general, the RMF model reproduces well the experimental data, and 
predicts a much subtle variance relative to other mass models. The neutron shell gaps, 
defined as A„(Z, A) = S2n{Z,A) — S2n{Z,A + 2), can be more clearly seen from Fig. [5l 
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which shows the shell gaps for N = 82 and 126 in the RMF approach in comparison with 
the data available and those in the FRDM, ETFSI-Q and HFB-13 models. The nuclei in 
the shadowed area are in the r-process path. At the = 82 shell, all the mass models 
except the FRDM model show a strong quenching effect (i.e., the shell gap drops) towards 
the neutron-rich side. The RMF shell gap is overestimated compared with the data available 
and it is around 2 MeV larger than other models for 45 < Z < 60. Regardless of this, the 
RMF model succeeds in predicting the enhanced double-magic effect at Z = 50 together 
with the HFB-13 model. For = 126, there is no sign of shell quenching observed in the 
r-process region for all the models. A unique feature of the RMF model is that it fully 
coincides with the available data and it is also the only model to reproduce the enhanced 
double-magic effect at Z = 82. In comparison, the other models fail to reproduce the trend 
of the known A^ = 126 shell. Towards the neutron drip line, the RMF prediction tends to 
enhance the shell until the maximum is reached around Z = QO while the other models have 
a roughly constant shell gap. 



III. SITE-INDEPENDENT r-PROCESS APPROACH 



Since the r-process is responsible for the synthesis of half the heavier nuclei beyond 



the iron group [sl, Q|, the basis of the nucleosynthesis mechanism have been extensively 
studied. Nevertheless, the location where it occurs has not been unambiguously identified. 
Current location candidates include the neutrino-driven wind off a proton-neutron star in 

neutron star mergers 



core collapse supernovae 



351, 
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40 



411 ]. jets in core 



collapse supernovae 4^, shocked surface layers of 0-Ne cores 43|], and gamma ray bursts 
{4^ . Because the specific astrophysical conditions among the different scenarios may change, 
solar r-process abundances 45] have been used in the past to constrain the astrophysical 
conditions using a site-independent approach [3, Q]- ^^i^ approach seed-nuclei (usually 
the iron group) are irradiated by neutron sources of high and continuous neutron densities 
Un ranging from 10^'^ to 10^^ cm~'^ over a timescale r in a high temperature environment 
(T ~ IGK). This superposition of r-process components (ra^,?") is needed to reproduce the 
overall shapes and positions of the solar r-process abundance. BaQa.d..e,u«e. 
to the exponential neutron exposures in the s-process 48|. The configuration of many r- 
process components seems to be also a reasonable approximation to the real r-process event. 
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For instance, one can think of it as the "onion" structure of neutron sources with different 
densities, where the seed-nuclei capture neutrons while moving through different zones with 
different thicknesses. In this paper we explore for the first time the application of the new 
RMF mass model to an r-process calculation and at the same time investigate the effect of 
nuclear physics uncertainty in the r-process. 

Due to the high neutron densities, neutron captures are much faster than the compet- 
ing /5-decays and an (n,7)<^(7,n) equilibrium is nicely established for every element. The 
abundance ratio of two isotopes in the timescale r can be expressed simply as 



where Y{Z, A) denotes the abundance of the nuclide (Z, A), Sn is the one-neutron separation 
energy, G{Z, A) is the partition function of nuclide {Z, A), and h, k and are the Planck 
constant, Boltzmann constant and atomic mass unit, respectively. Neglecting the difference 
in the ratios of the partition functions and the atomic mass, one can easily see that the 
isotopic abundance distribution P{Z, A) and the abundance maxima in each isotopic chain 
are determined by n„, T and Sn- Approximating Y{Z,A + 1)/Y{Z,A) ~ 1 at the highest 
isotopic abundance for each element, and all other quantities being constant, the average 
neutron-separation energy Sn, calculated by 



is the same for all the nuclides with the highest abundance in each isotopic chain. Tg 
denotes the temperature in 10^ K. Higher temperature or lower neutron density will drive 
the r-process path towards the valley of stability. Due to the pairing correlation the most 
abundant isotope always has an even neutron number A^. 

If fission is neglected, the abundance flow from one isotopic chain to the next is governed 
by /9-decays and can be expressed by a set of differential equations: 




exp 



Sn{Z,A+l) 



(3) 






A 




(5) 



A 
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where X^'^ is the total decay rate of the nuchde (Z, A) via the /3-decay and the delayed 
neutron emission, and Y{Z) = "^^Y^Z, A) = "^^P^Z, A)Y{Z) is the total abundance 
in each isotopic chain. Using Eqs. ([3]) and ([5]), the abundance for each isotope can be 
calculated. After the neutrons freeze out, all the isotopes will proceed to the corresponding 
stable isotopes via /3-decays. 



IV. CALCULATIONS 



In the present calculation, unknown one-neutron separation energies S'„ were calculated 



from the RMF approach 
experimental data 
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3l| and /5-decay properties were taken from Ref. 49|]. Available 



50| was used when available. Similar to the method used in Refs. 
46|, |47[, we applied sixteen components with neutron densities in the range of 10^''-3 x 10^^ 
cm~^ in our calculation. We chose a temperature T = 1.5 GK. We assumed that for this 
temperature the irradiation time r and the corresponding weight w follow the exponential 
dependent of neutron density n„, i.e.. 



w;(n„) = n", r(n„) = bxn 



c 

n ' 



(6) 



where a, b, c are parameters to be fixed. These parameters can be obtained from a least- 
square fit to the solar r-process abundances. We further assume that the longest neutron 
irradiation time has to be longer than 0.5 s but shorter than 20 s. The exponential relations 



in Eq. ([6]) have been observec 
and chronometers studies 



when fitting the three r-process peaks 



47 



52|. 



5l[ | and used for stellar 



It was found that r-process components with r(?7,„) = 0.454n^'°^° s and w{nn) = 2.1n^^'^ 
best reproduce the solar r-process abundance. Fig. shows the contribution of the four- 
group weighted r-process components after /9-decays to the resulting best fit. The black 
solid curve with isotopic abundances normalized to A = 130, corresponds to the fit using 
all the sixteen components. The green, red, blue and grey dashed curves are the sum of 
the abundances calculated with log(n„) ranging from 20 to 22.5, 23 to 24,5, 25 to 26.5 and 
27 to 27.5, respectively. The first six components with log(r2„) between 20 and 22.5 seem 
to account for the A = 80 abundance peak. The four components with log(?7,„) between 23 
and 24.5 are responsible for the overall structure of the r-abundance curve beyond A = 120 
and the remaining components only improve the description of the theoretical calculation 
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for A > 150. In general, the fit is found to reproduce well the solar r-process abundances 
and also the position of the abundance peaks. 

The r-process runs relatively close to the /9-stability line around the shell closure, thus the 
experimentally known mass values around = 82 shell significantly influence the abundance 
distribution after the second abundance peak. Taking our best simulation using the RMF 
masses as an example, the ratio between the abundance at A = 130 and the abundance at 
A = 195 is 2.8 when taking the experimental data into account, and it increases to 24.6 
if the experimental data are not used. The same ratio decreases from 4.0 to 1.3 for the 
best simulation using the FRDM masses (to be discussed below). In order to minimize the 
contribution from the theoretical uncertainty of known masses, experimental information 
was included in the calculations. 

In order to investigate the impact of theoretical uncertainty of unknown masses in an 
r-process calculation, we also performed the same procedure using instead the FRDM, 
ETFSI-Q and IIFB-13 mass predictions while keeping the same /3-decay properties. The 
astrophysical conditions determined by using various mass inputs are shown in Fig. [71 The 
obtained superpositions of sixteen r-process components for all the mass tables are collected 
in table [B Similar to the astrophysical conditions obtained from ETFSI-Q simulations, the 
astrophysical condition using the RMF mass input requires a relatively constant weight- 
ing factor for different neutron densities. However, the FRDM and HFB-13 cases favor a 
large weighting factor for the low neutron density. As for the neutron irradiation time, the 
best RMF fit requires component durations of as long as 6 s while the FRDM and HFB-13 
simulations only requires up to 1.5 s. The ETFSI-Q component durations are somewhat 
in between. Moreover, it may be worth mentioning that the simulations using the FRDM 
and HFB-13 masses demand almost identical astrophysical conditions. Using a lower tem- 
perature T = 1.35 GK, a similar calculation based on the ETFSI-Q mass model is carried 
out in Ref. [47 1. As shown in Fig. [3, their obtained neutron irradiation times are in good 
agreement with our calculation using the ETFSI-Q masses and FRDM half-live inputs, but 
the weighting factors differ. The superposition obtained in Ref. 47| demands a more sharp 
evolution of the weighting factor as a function of n„. Since a lower temperature of T = 1.35 
GK in our calculations only weakly impacts the condition obtained, the difference should be 
due to the different /5-decay properties used in that work. Based on a full dynamical network 
calculation, faster time-scales of the order of hundreds of milliseconds, are found in Ref. 53 1 
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for an r-process in the neutron-wind scenario of core-collapse type II supernovae. However, 
this different time scale can be at least partially attributed to different seed nuclides. In 
their calculation, the r-process stars from a seed distribution containing neutron-rich nuclei 
with mass numbers between 80 and 100, while ours starts from ^^Fe. 

Calculated solar r-process abundances after /5-decays using different mass models are 
displayed in Fig. [H Shadowed areas show the regions with underproduced abundances before 
the neutrons freeze-out. After /^-decays to the stability line, these gaps are too large to be 
completely filled in by /^-delayed neutron emissions. It should be pointed out that the solar 
system r-process abundances are defined as the abundances not produced in the s-process 
and p-process that still have to be created elsewhere to explain the solar system abundances. 
Although it is thought that the r-process is responsible for the majority of those isotopic 



abundances with Z > 56, its contribution to the lighter elements is still debatable |57l . |58 |. 
It is possible that some of the discrepancies in the reproduction of the low mass abundances 
may be due to an additional nucleosynthesis component creating some of those abundances. 
However, since astrophysical conditions and nuclear properties both affect the resulting r- 
process abundances, one need to determine or at least understand the uncertainty in the 



591 . In this paper 



nuclear physics properties in the future works to disentangle both effects 
we only discuss possible nuclear physics reasons for such underabundances. 

-R-process abundances calculated with all nuclear mass models result in abundance un- 
derproduction at A ~ 120 and A ~ 170. Traditionally, the underestimation of the isotopic 
abundances before A ~ 130 peak has been attributed to the overestimated strength of the 



= 82 shell closure 



471, 



511 ] in the theoretical nuclear physics model even though the 



experimental evidence is still debated [5J, l55|, l56| . Since it is not possible to do a complete 



study of the shell-quenching effect at the single particle level which should affect more nuclei 
than the one with A^ = 82, we only study the effect of a reduced shell closure by artificially 
decreasing the shell gap energies at A^ = 82 in the RMF and FRDM models by 2 MeV and 
1 MeV, respectively. In such a way, the shell gaps interested for the r-process would roughly 
have the same values as those in the quenched models ETFSI-Q and HFB-13. Eventually, 
a better agreement with the observation at A ^ 120 is obtained as shown in Fig. [8]|^a-b). 
This can be easily understood as follows. A reduction of shell gap leads to a nuclear matter 
repopulation in the isotopic chain according to Eq. [31 i?-process waiting points located at 
A^ = 82 move closer to the valley of stability and thus some of the underabundance can be 
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filled. Furthermore, based on Fig. [3] one could expect that the quenched shell gap at = 82 
would not affect the abundance around A ~ 115. It is interesting to note that the r-process 
simulation using the shell- quenched ETFSI-Q model in Fig. M^c) show a good agreement 
with solar abundances pattern at A ~ 120 together with a large underproduction at a lower 
mass number A ~ 115. 

The abundance trough at A ~ 115 for the FRDM, RMF and ETFSI-Q models can be 
related to the additional bump oi Sn a-t A = 110 — 120 in Fig. [3l and thus the associated 
nuclear shape transition. In the case of the ETFSI-Q model, nuclear shape changes for 
prolate to oblate and then to spherical nuclei with = 82. This transition leads to a 
deviation from the approximate relationship between neutron separation energies and mass 
number for each isotope, and can be clearly recognized in Fig. [31(c) by the sudden increase 
of the separation energies. In order to see the sensitivity of the r-process calculation to the 
effect of the nuclear shape transition, we lowered the separation energies of ^^^'^^°Mo by 1 
MeV in the ETFSI-Q model, but kept the other nuclear physics input unchanged. Those 
isotopes are in the r-process path and show a bump in the one-neutron separation energies as 
a function of mass number. As shown in Fig. [8](c) the A ~ 115 trough is largely filled in. A 
similar analysis shows the same conclusion for the FRDM and RMF mass models. Although 
we have mainly focused in the underabundance below A^ = 82, similar conclusion can be 
drawn for the trough around A ~ 170. As an example. Fig. ^c) shows that the trough is 
almost completely filled in by lowering the separation energies of ^^^Pm and ^^^Sm by 1 MeV. 
This suggests that the potentially wrongly assigned location of the shape transition before 
the neutron magic number in the theoretical predictions can lead to the troughs before the 
abundance peaks. 

Of all the mass models, the HFB-13 is the only one that shows a smooth one-neutron 
separation energy change from Sr to Ru. As a result, the r-process waiting points are contin- 
uous and there is not apparent gap in the r-process path (see Fig. [3]^d)). Only modifications 
in the nuclear masses would not result in the filling of the the abundance gap at A ~ 115. 
Such underproduction may be traced back to the /3-decay properties. By increasing the 
/3-decay half-lives of the critical nuclei ii3,ii5,ii7y Iqj gye times, we found that the trough 
before the A ~ 130 peak in the HFB-13 case can be nearly filled in as shown in Fig. [8](d). 
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V. SUMMARY 



We have applied the most recent comprehensive mass models, the non-relativistic 
microscopic-rooted HFB-13 and the relativistic RMF in r-process calculations. For the 
sake of comparison, we also included the widely used macro-microscopic models FRDM and 
ETFSI-Q. Of these models, the HFB-13 and RMF models are used for the first time in such 
calculations. Based on a simple r-process model, it is found that all mass models reproduce 
the main features of the solar r-process pattern and the position of the abundance peaks. 
Since r-process simulations have to rely on predicted nuclear physics properties of unknown 
regions in the nuclear chart, we have compared the predictions of different mass models. We 
have also made a systematic study of the influence due to the mass model uncertainty in the 
application of the r-process and thus in the required astrophysical conditions. This nuclear 
physical uncertainty is very important for the complete understanding of the r-process since 
the results of more modern full dynamic r-process calculations depend on the nuclear mass 
input used. It is found that the deduced astrophysical conditions like the neutron irradia- 
tion time of the r-process can be significantly different depending on the mass model used. 
Among the different models, the simulation using the RMF masses requires a longer time 
scale (up to a factor of 4) than those using FRDM and HFB-13 models. Furthermore, it is 
found that the optimal astrophysical conditions obtained using the ETFSI-Q and RMF mass 
models require a relatively constant weighting factor for neutron densities in the range 10^^ 
to 10^^ cm~^, while the FRDM and HFB-13 simulations favor a large weighting factor at 
low densities. In addition, we have explored the possible deficiencies in different mass mod- 
els, and found that the observed abundance underproduction before the abundance peaks 
in all the models can be a combined and complex effect of both shell structure and shape 
transition. An exception is the underproduction at A ~ 115 in the HFB-13 model which 
can be attributed to incorrect /5-decay rates. Future experiments are needed to determine 
the strength of the shell closure towards the neutron drip line as well as the precise locations 
of the shape transition toward the shell-closures. 
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FIG. 1: (Color online) The rms deviation arms of one-neutron separation energy Sn with respect 



to experimental data 



33l | as a function of the distant from the /3-stability line e = Zq — Z for 



different mass models, where Zq stands for the proton number of the most stable isotope in the 
isobaric chain with mass number A. 
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FIG. 2: (Color online) The differences between one neutron separation energies S„ predicted in 
the RMF model and those in the HFB-13 model. The magic proton and neutron numbers are 
indicated by pairs of parallel lines, and also the present border of the data with known masses are 
shown by solid lines. 
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FIG. 3: (Color online) The average one-neutron separation energies around the N = 82 shell in the 
RMF model, in comparison with those in the FRDM, ETFSI-Q and HFB-13 models as a function 
of mass number A. For simplicity only nuclei with even N are plotted. The corresponding r-process 
paths calculated using different mass inputs are also indicated by dots, and labeled here are those 
isotopes with more than 10% population of each isotopic chain. 
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FIG. 4: (Color online) Same to Fig. [3] but around the N = 126 shell. 
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FIG. 5: The neutron shell gaps A„(Z, A) = S2n{Z, A) - S2n{Z, A + 2) for = 82 and 126 in the 
RMF approach compared with those in the FRDM, ETFSI-Q and HFB-13 models together with 
the data available. The nuclei in the shadowed areas are involved in the r-process paths based on 
our calculations. 
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FIG. 6: (Color online) The effect of various weighted r-process components on the resulting fit 
after /^-decays in the best superposition using the RMF masses. The calculated total isotopic 
abundances (in the logarithm scale)) is normalized to A = 130. 
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FIG. 7: The best configuration of sixteen r-process components which reproduce the solar system 
r-process abundances with different mass inputs at the temperature T = 1.5 GK. The neutron 
density is in the unit of cm~^. The weighting factor uj{nn) and the neutron irradiation time 
T{nn) are shown in the upper and lower panels as a function of neutron density n„. In the upper 
panel, the weighting factors for the FRDM are completely overlaid by those for the HFB-13 models. 
The fit from Ref. [47|] is also plotted for comparison. The total weighting factor has been normalized 
to 100. 



TABLE I: Our best fits to the solar r-process abundances for different sets of nuclear mass models. 
The first column is the mass model employed. The last two columns are the weight uj and the 
relevant neutron irradiation time r (in unit of second), respectively. 

Mass model uj t [s] 



RMF 


2.1xn0-020 


0.454 xnO-040 


FRDM 


3.0E4xn-0-i6i 


0.013 xnO-O'^s 


ETFSI-Q 


54.4xn-0-040 


0.499 xn°-°25 


HFB13 


2.8E4xn-0-i60 


0.007 xnO-085 
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FIG. 8: (Color online) Our best fits to the solar r-process abundances (in the logarithm scale) 
under different nuclear mass inputs. The /3-decay properties are taken from the FRDM model (49l |. 
The best fits are displayed as blue-solid lines. In the sub-figure (a), the red-dashed curve is the 
same as the blue-solid curve but with a shell closure at = 82 2 MeV smaller. In the sub-figure 
(b), the red-dashed curve is the same as the blue-solid curve but with a shell closure at = 126 
1 MeV smaller and separation energies of ^^^Pm and ^^^Sm 1 MeV smaller. In the sub-figure (c), 
the red-dashed curve is the same as the blue-solid curve but with separation energies of ^^^Pm and 
ii8,i20gj^ 1 MeV smaller. In the sub-figure (d), the red-dashed curve is the same as the blue-solid 
curve but with half-live of isotopes ii3,ii5,ii7y five times larger . The shadowed areas correspond 
to the range where the abundances of these isotopes are largely underestimated before neutrons 
freeze out. 
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